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The exact solubility of the Kitaev-type spin honeycomb lattice model was proved by means of a 
Majorana fermion representation or a Jordan- Wigner transformation while the explicit form of the 
anyon in terms of Pauli matrices became not transparent. The nonabelian statistics of anyons and 
the fusion rules can only be expressed in indirect ways to Pauli matrices. We convert the ground 
state and anyonic excitations back to the forms of Pauli matrices and explicitly demonstrate the 
nonabelian anyonic statistics as well as the fusion rules. These results may instruct the experimental 
realization of the nonabelian anyons. We suggest a proof-in-principle experiment to verify the 
existence of the nonabelian anyons in nature. 
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Introduction The Kitaev-type spin honeycomb lattice 
models have attracted many research interests for the 
possible nonabelian anyonic excitations in these exactly 
soluble two-dimensional models [l| . There are two topo- 
logical phases for these kinds of models. The topologi- 
cally trivial A phase is an abelian anyon phase which is 
equivalent to that in Kitaev's toric code model The 
topologically nontrivial B phase is within the same uni- 
versality class of the Moore-Read Pfaffian states in the 
fractional quantum Hall state Isl kI and the vortex exci- 
tations are nonabelian anyons [sl, H, 0] . 

In solving these kinds of models, a key technique is the 
usage of the Majorana fermion representation of the spin- 
1/2 operators, either via Kitaev's Majorana fermions or 
the Jordan- Winger transformation However, the 

shortcoming to introduce these Majorana fermions is that 
the ground state and the elementary excitations are hard 
to be expressed by the original spin operators, i. e., Pauli 
matrices. Then the nonabelian fusion rules and statistics 
may not be directly shown in Pauli matrices' language 
[lO|. Meanwhile, experimentally exciting, manipulating 
and detecting anyons may be more practical by using 
the spin operators, as recently s ugg ested or done for the 

and for the A 
Therefore, to 



toric code model [ll|, [13, 0, lU, 
phase of Kitaev honeycomb model [l^, |l£ 
explicitly demonstrate the nonabelian anyonic statistics, 
one needs to express the ground state and elementary 
excitations in the spin operators. Chen and Nussinov [9| 
have studied a real space form of the ground state of the 
Kitaev honeycomb model and applied it to the A phase 
with abelian anyons. However, for the more interesting 
B phase, it was not figured out yet. 

In a recent work, one of us, with Wang, has provided 
a generalized model of the Kitaev honeycomb model by 
adding three- and four-spin couplings [j, [20| . With this 
generalized Kitaev-type model, we showed the equiva- 
lence between the B phase with the breaking of the time 
reversal symmetry and the Moore- Read Pfaffian state. In 
that work, we map the model in the honeycomb lattice 



to a spinless fermion model in a square lattice. The task 
of the present paper is mapping back the ground state 
and anyonic(vortex) excitations obtained in the square 
lattice to their honeycomb lattice version. 

After writing down these states in spin operators, i.e., 
Pauli matrices, we can check the nonabelian statistics of 
anyons and fusion rules of the excitations. As expected, 
these results agree with those in the previous abstrac- 
tive study in Kitaev's original work [ij. Moreover, these 
explicit forms with well-known Pauli matrices may help 
the readers who are not in this special field to understand 
those formal descriptions made by Kitaev. It may also 
instruct the experimentalists to realize these states and 
verify the existence of the nonabelian anyons in nature. 
Similar to the case of the abelian anyons in the toric code 
model 14, 3, l3], we may also expect to excite, manip- 
ulate and detect the nonabelian anyons in atomic spin 
systems in optical lattice. We may design a proof-in- 
principle experiment to show the existence of nonabelian 
anyon in nature by means of the techniques proposed and 
developed recently in the photon graph state [lll,[ll[ll, 
and a nuclear magnetic resonance system [l5| . The min- 
imal lattice needs only six sites, which is accessible to 
the current experiments, as several experimental groups 
have done to the toric code model [ll[l3,[ii. 
Model and ground state We consider a Kitaev-type spin 
model in a honeycomb lattice with a three- and four-spin 
couplings 0, 1131 
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where cr^'^'^ are Pauli matrices, a;-,j/-,z-links are shown 
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in Fig. [Tl 'w' and 'h' label the white and black sites of 
lattice, and ex^Sy, are the positive unit vectors, which 
are defined as, e.g., 



a2 1 ^aia^ 



Jx,y,zy >^ and Xj. y are tunable real parameters. This is a 
system with the breaking of the time-reversal symmetry. 
There is a Z2 gauge symmetry generated by 

W^P<.=««« (2) 
where Pa labels a plaquette (see Fig. [1]). That is. 



According to Lieb's theorem [2l|, the ground state of 
the system is Z2 vortex free state. This means that 
ibbbb+e^ = 1 in the ground state sector. Following the 
track in refs. finally, the Hamiltonian in the ground 

state sector is given by 

Ho = Jz J2^dld, - 1/2) + ^KJ2id^d^+e.^ + 44+eJ 
i i 

+Jx{d\di+e^ - did\^^J + Jy{d\d,+ey - d^dl^^J (5) 

+ Xx - di+eji) + ><v^{d\+eyd\ - di+e^di) 




FIG. 1: The honeycomb lattice and plaquette. The thick links 
define Aj and also the proof-in-principle configurations. 

We first recall how to work out the ground state of 
this model. The first step is introducing the Majorana 
fcrmion representation of the spin operators cr" Jl| or em- 
ploying the Jordan- Wigner transformation 0, d, 3| ■ For 
convenience, we use the Jordan- Wigner transformation 
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where i lables a z-link and tpi^ and are Majorana 
fermions which satisfy "i/ji^ i/jj-j — —^j^ipi^, bi^bj^ ~ —bj^bi^ 
for is ^ it] and ij^i^ ~ ^^jt ~ ^^'^ ^'^^ anti- 

commutative. The order of the sites is defined as follows: 
is > it if the zig-zag horizontal line that is belongs to is 
higher than that of jt or if is is on the right hand of jt 
when they are in the same line. 

After this Jordan- Wigner transformation, the Hamil- 
tonian reads 
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(4) 
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where Jx^y = 



and Xx^y = 



The spin- 
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less fermion di = ■^{ipi^ + is located at the z-link 
and all z-links form a square lattice. Taking A^; -I- = 
and Ay — k, this Hamiltonian describes a, Px + ipy-wstve 
pairing state in this square lattice. As we have shown, 
the phase diagram consists of two phases: the topologi- 
cally trial A phase and non-trivial B phase. The phase 
boundary is the lines: J2 i JxT Jy if we restrict to > 
and Jx,y > 0. The A phase is an abelian anyon phase 
which is equivalent to that in the toric code model, which 
was studied before [H, Q . If we do not consider the high 
energy Majorana fermion excitation, the effective Hamil- 
tonian of the A phase in the honeycomb lattice reads [H 



= -J, 



(6) 



z — links 



where Aj — ipiip^tpQip^ (see Fig. 1) with Xx^y = k = 

and Jeff = ^Qj3 for j, <C Jz- One may directly check 
that the ground state is given by 



\GA)^Y[(.l + A,)Y[il + WpM), 



(7) 



because Wp\Ga) = Aj\Ga) = \Ga) with |(/,) = | t • ■ • t) 
a reference state. Here the Majorana fermion excitations 
have a high energy 2Jz and has been neglected. This 
Hamiltonian has a Z2 x Z2 gauge symmetry generated 
by Wp and Aj with A^j ^ ^ 1 and [^j,Wp] = 
0. We know that Aj and Wp are corresponding to the 
'electric charge' and 'magnetic charge' in the toric code 
model. The low energy excitations have Aj — —1 or 
Wp = — 1, which are e and m vortices in the toric code 
model and obey the mutual semion statistics [Ij]. It was 
noted that the fermion excitations may not be ignored in 
some braiding processes |l9l| and there is a controversy 
to this matter recently |22| . 

In the B phase, we do not have a conserved 'electric 
charge' Aj. Since the ground state sector in the continu- 
ous limit is a Px + ipy-waye BCS theory, the ground state 
in the square lattice may be written down, which is 0] 



\Gb) = ( 



1 



-44)^/2 = pf(- 
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where Zj = jx+ijy is a complex number with j = {jx, jy) 
the lattice site label. The vacuum state |0) is defined 
by dj\0) = while the state \D) = Jli^lO) satisfies 
dj\D) = because the square lattice is filled. The vac- 
uum state has been written back in terms of Pauli matri- 
ces 0, which is given by |0) oc llji'^-Q3)Ylpi^+WpM) 
due to dj{l — Qj) = where Qj = iipj^ipji^. 1 + Wp fac- 
tor is introduced because Wp{l + Wp) = 1-1- Wp ensures 
the vortex free of the ground state, i.e., Wp\Gb) = \Gb)- 
Because oi (^^{l + Qj) = 0, djiZ?) =0, which is equivalent 
to Qj\D) = ID). Therefore, 



\D)^\{{l^Q,)\{{l^Wpm. 



(9) 



Excitations A Majorana fermion excitation on the ground 



state is given by ^j^\GA,B)- Due to [Wp,^jJ 



0, 



ipjA^A.s) is also a vortex free state. According to the 
original Hamiltonian, the energy cost to excite a Majo- 
rana fermion is given by 

E^-Eg = {GAMi>js\HHjs\GA,B) -Eg = 2J„ (10) 

where Eg is the ground state energy. Note that ipi^ relates 
to by ^l^i^lGs) = -ii>i„\GB)- 

We are now going to create the vortex excitations 
which are defined by Wpap\GA,B) = —<^p\Ga,b) and 
Wpi(Jp\Ga,b) = ctpIGa.b) for P' ^ P. Two operators 
obey these requirements: 



.(1) 



.(11) 



The vortex is located at the plaquette P with if, being 
its 'ai' (See Fig. 1). The sites i;, — 1, if, — 2, • ■ • are also 
marked in Fig. 1. One may also define a vortex at the 



same plaquette through a 



p 



• • • and 



err o- 



Creating a single vortex costs an infinite energy, e.g., 
in the B phase 

{GB\of\H\G'^^\GB) ^Eg + 2(J, + Jy){ih - (-oo)Xl2) 

Therefore, it is impossible to excite a single vortex. How- 
ever, exciting a pair of vortices spends finite energy which 
is dependent on the difference of the site labels of two 
vortices, e.g., two adjacent vortices (crf^CTj? _2crf^_4 • • •) • 
(<+2«-2---) = <-H2> which costs energy 4(J^ -I- 
Jy). In the A phase, the energy cost of a pair of cr'^^^ is 
2Je//, which is much lower than the energy cost to ex- 

(1) (2) (2) (2) 

cite ip as Jz ^ Jx,y Exciting pairs Cp Vp/ and Cp Vp/ 
costs energy 2(2Je// -t- Jz) and ^{Jeff + Jz), which are 
also the high energy excitations. 

Fusion rules To see nonabelain anyonic fusion rules, we 
focus on the B phase to study the fusion rules. Since 
ipf = 1, the fusion rule of the Majorana fermions is ip x 
— 1. For the vortex excitations, one has 



4i).4i) = l, 42). ^(2)^^^ 4i).a(?)oc^, 

; (1) (2) ; (2) (1) 



Define two vortex operators ap = aap"^ +Pcrp^ and dp ~ 

a*cr^^^-K/3*cr^^\ which obey apa], = \a\^ + \p\^ + i{af3* - 
a*(3)ipi^. We can not distinguish vortex pairs apapi, 
apffp,, ap(jp, and CpCTp' because they are energetically 
degenerate. This means the equivalence between ap and 
(Tp and the fusion rule for the vortices isaxa = l + ip. 
On the other hand, the paired vortices are described by 
a Pfaffian wave function Q] 



Pf 



{Zj - Z,p){Zj, - z, ) + {j ^ j') 



(14) 



no matter what kind two vortices are. This also implies 
the equivalence between ap and a p. Summarily, the fu- 
sion rules in the B phase are 



■0x?/' = l, a X a — 1 + Tp, ^ X a 



(15) 



These nonabelian fusion rules are the same as those in 
the Ising model. All fusions cost energy in the same order 
as that to create a ip and a pair of vortices. 
Braiding matrix The Majorana fermions are anti- 
commutative which gives Rf^ = — 1. To see the braiding 
matrix elements for the vortices, we rotate vortices coun- 
terclockwise. In the B phase, we consider two vortices at 
P2 and P3 in Fig. 1. We have two ways to rotate them. 
One way is moving the vortex at P2 to Pi first, then mov- 
ing the vortex at P3 to P2, and moving that at Pi to P3 
(See Fig. 2(up panel)). This exchanges two vortices and 
the three steps is given by acting crj^ • a 



af m turn 

Jo 



on ap^ap^lGB)- A single ct^'^ action means, while one 
vortex is moving, (say ap^ — > ap-^) another one, say crpg, 
becomes (Tp =i})xap^. Thus, this exchange is accompa- 
nied by non-trivial fusions or creation and annihilation 
of the Majorana fermions and then the braiding matrix 
element is denoted by R?','^ = a^ a^ a'^ = i. Rotating 

Jb jh Jb 

vortices clockwise, we have i?^""^ = '^jb'^jt'^jb ~ 
This also means R'^'f' = R^" = -i. Another way 
to exchange two vortices are moving the vortex at P2 
to Pi and that at P3 to P4 simultaneously, then mov- 
ing that at Pi to P3 and P4 to Pi in the same time 
(See Fig. 2 (low panel)). This exchange corresponds 
*o ^i^^Jb^ji'^l'^P^o-ftlG's) = f^tcr|,CTp2crp3|Gp) = 
(TpjCrpglGs) since flt^l^ IGp) = \Gb)- Simultaneous 
(tJ^ and (tJ^ action creates two Majorana fermions ipj^ 
and ipj^ but the relation iijji,..ijji^ IGpi) = \Gb) means no 
fermion is created at any stage [l9| . This defines a braid- 
ing matrix element R.i°' = 1 to exchange. Therefore, the 
nonabelian braiding matrix for the B phase is given by 

Rf'^' = _1, i?^'^ = 1, ^ jif<^ = jio^<^ = j. (16) 

Here we do not put in the abclian phase factor e~"^/® 
which comes from the Pfaffian (fT4|) . Restoring this phase 
factor, one has R^'' = g-'-^/s and R'^" = e'^''/^ and the 
braid matrix are the same as that of the Ising model 



(13) [ESS- 
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FIG. 2: Exchange of the vortices. ' — ' at a plaquette labels a 
vortex. Pauli matrix above the filled arrow means the hermi- 
tian rotation operator acting on the white site '5' while below 
the filled arrow means the operator acting on the black site 
'2'. Up panel: The exchange leads to a phase factor i for the 
B phase. Low panel: a bosonic exchange for the B phase. In 
the A phase, the up panel is forbidden . 



Conclusions We have converted the ground state and 
elementary excitations from the fermionic representation 
in a square lattice to the original spin representation in 
the honeycomb lattice for the Kitaev-type model. Pauli 
matrix version of these states leads to an explicit demon- 
stration to the non-abelian statistics of anyons in this 
model. We showed the nonabelain fusion rules and cal- 
culated the non-abelian braiding matrices. We proposed 
a proof-in-principle experiment to create, manipulate and 
detect the nonabelain anyons in nature. 

The authors thank Ville Lahtinen and Julien Vidal for 
the helpful comments. This work was supported in part 
by the national natural science foundation of China, the 
national program for basic research of MOST of China 
and a fund from CAS. 



Experimental implications Recently, there are several 
proposals to excite, operate and observe the abelian 
anyons in the toric code model. Most of them are based 
on the atoms or molecules in optical lattice 0,111 [13, 
Sin. Since the nonabelian braiding matrix (|16p is not 
related to the Pfaffian factors in eqs. ([5]) and all 
the states involved in do not have the site-dependent co- 
efficients if we neglect the Pfaffians. Thus, all techniques 
applied to the toric code model may be employed to the 
present model. For example, load cold atoms in a hon- 
eycomb optical lattice and manipulate an ancillary atom 
as proposed in Ref. 16]. 

We can also design a possible proof-in-principle exper- 
iment to the nonabelian anyons by using the systems 
to prove the abelian anyons in the toric code model 
11, [13, [13, 15 1 . As Han et al designed a scheme to 



demonstrate the abelian anyons in the toric code model, 
the minimal lattice needed to rotate or exchange two vor- 
tices are six sites connected, e.g., by the the thick links 
in Fig. 1. For this minimal lattice, we prepare the state 

1^) 



\D) cx (1 -|- ai^aZa^a 
(l + aK)(l 



4"5 



(17) 



This is an entangled state of 28 pure states and a bit 
complicated to be prepared experimentally but is still 
accessible. The two vortices state may be given by, 
e.g., crllD) which creates two vortices at P2 and P3 
because {cf, erf frffTlcrg } = {o"!, erf o-|cr|crf } = and 
[cr|, o-fcrlcrf ] = 0. Two exchanges described before, 
erf (T|crf fjf and cr2'^2^2j result in the signs ±, respectively. 
A full nonabelian two vortices includes a Pfafhan factor 
(fT4)l . which contributes an abelian phase factor. We do 
not know if it is possible or not to prepare such a state. 
Fortunately, to see the nonabelian braiding matrix p6p . 
one needs to simply act a two vortex operator on \D) and 
check those two different exchanges are enough. 
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